Proceedings of the DPF-2011 Conference, Providence, RI, August 8-13, 2011 



1 



Search for Flavor Changing Non-standard Interactions with the MINOS 

Experiment 



O 

(N 
-i— > 
O 

O 
o 



X 

in: 

D 
43 



> 

o 
o 

<3\ 



Zeynep Isvan, for the MINOS Collaboration 

Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA, USA 

MINOS searches for neutrino oscillations using the disappearance of muon neutrinos between two 
detectors, over a baseline of 735 km. We recently reported the most precise measurement of neutrino 
oscillations in the atmospheric sector and the first tagged measurement of antineutrino oscillations. 
The neutrino mass splitting and mixing angle are measured to be |Am 2 | = 2.32t ( Q'J 2 , 
and sin 2 26 > 0.90 (90% C.L.) for an exposure of 7.25 x I0 20 protons-on-target (PoT). Antineu- 
trino oscillation parameters are measured as Am 2 = (3.36ig'.4o(stat.) ± 0.06(syst.)) x 10~ 3 eV 2 and 
sin 2 (2#) = 0.86to:i2( stat O ± 0.01(syst.) with an exposure of 1.7 x 10 20 PoT in NuMI antineutrino 
running mode. We use the apparent difference in neutrino and antineutrino oscillation parameters 
to constrain non-standard matter interactions which could occur during propagation through the 
Earth's crust to the Far Detector. 



I. INTRODUCTION 

MINOS observes a difference between neutrino and antineutrino standard oscillation parameters Q], Hj]. The 
statistical probability that the neutrino and antineutrino data have the same oscillation parameters is 2% 0, Q . 
So-called non-standard matter interactions (NSI) of muon neutrinos and antineutrinos with the Earth's crust 
during their propagation can be a possible explanation of this difference, since neutrino (matter) and antineutrino 
(antimatter) survival probabilities are altered in opposite directions by this effect. MINOS is particularly well 
suited to perform this measurement since its detectors are magnetized and neutrinos and antineutrinos can be 
distinguished on an event- by- event basis. 

This document details the non-standard interactions (NSI) analysis which uses a 7.09 x 10 20 POT neutrino 
sample and a 1.71 x 10 20 POT antineutrino sample. In this analysis neutrino and antineutrino vacuum oscillation 
parameters, Am 2 and sin 2 29 are assumed identical and the entire difference between the observed rate of 
disappearance is attributed to non-standard interactions. 



Neutrino Oscillations 



Neutrino oscillations occur because neutrino flavor eigenstates are not identical to the mass eigenstates but 
instead are related by the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix, U. The probability that 
a neutrino of energy E created with flavor a will be detected in a flavor state j3 after traveling a distance L is 
given by 
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where rrij is the mass of the jth mass eigenstate. The final expression for oscillation probability depends not on 
the absolute masses of the neutrinos but on the squared difference between masses, Arn 2 2 and Am 2 ^. Because of 
the large difference between the magnitudues of these two mass squared splittings, experiments such as MINOS 
are sensitive to mixing between two of the three flavors. In the two-flavor approximation the survival probability 
for a muon neutrino is given by: 



P{v u -+ Vv ) = 1 - sin 2 2(# 23 ) sin 2 (l .27^f^j (2) 

This flavor change is a quantum mechanical phenomenon observed in vacuum. In the presence of matter along 
the neutrino's path, this probability is altered since neutrinos interact with the matter and scatter. While all 
flavors scatter via neutral current interactions, this effect is larger for electron neutrinos because they can also 
partake in charged current interactions with eletrons. These standard matter interactions affect the survival 
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probability of the electron neutrino to a greater extent than that of muon and tau neutrinos, since muon and tau 
particle densities in normal matter are miniscule compared to electrons. There may be, however, non-standard 
matter interactions that influence survival probabilities of muon neutrinos. 

B. Non-standard Interactions 

If neutrinos participate in non-standard interactions, the flavor Hamiltonian will receive contributions from 
this effect similar to the standard matter interactions that give rise to the Mikheyev-Smirnov-Wolfenstein 
(MSW) effect 0,0. In general, NSI can be flavor changing or flavor conserving, with amplitudes proportional 
to the standard MSW matter effect £nnV, e TT V, and en T V, with V = ^/2GpN e . 
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MINOS is most sensitive to muon neutrino disappearance and the vacuum oscillation measurements assume 
a two-flavor model with transitions into tau flavor neutrinos. The limits on are already stringent and the 
difference between 77 M and oscillations is more sensitive to e Mr than e rT therefore the flavor-conserving 
part is neglected in the following discussion. Only the real part of e^ T distinguishes neutrinos from antineutrinos; 
the CP-violating imaginary part is neglected. This gives the following two-flavor NSI Hamiltonian, where e^ T 
is the real-valued, flavor changing neutral current contribution to non-standard interactions: 
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Using the fact that V —> —V for antineutrinos we can write the full Hamiltonian, H — Hq + H matte r as: 



1 sin 23 cos 23 ± e^V COS 2 23 



(5) 



with which the survival probability becomes, in units of GeV for neutrino energy, km for baseline and eV 2 for 
mass splitting: 
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Figure [T] shows a comparison of the vacuum survival probability to the probability with non-standard matter 
interactions at maximal mixing. 
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FIG. 1: Survival probability of neutrinos and antineutrinos in the presence of non-standard matter interactions upon 
traveling 735km from the source. MINOS best fit values for standard oscillations are assumed, Am 2 = 2.32 x 10~ 3 eV 2 
and sin 2 26 = 1.0. 



MINOS is a long baseline neutrino experiment which searches for neutrino oscillations primarily in the atmo- 
spheric sector by sending a beam of predominanly muon neutrinos 735 km from Fermilab to the Soudan mine 
in Minnesota. Two functionally identical detectors measure the neutrino energy spectrum before and after os- 
cillations occur. This allows for a meeasurement of neutrino oscillations by looking for a disappearance between 
the Near Detector, 1km downstream of the target at Fermilab and the Far Detector, 735 km downstream at 
the Soudan Mine. 

The NuMI (Neutrinos at the Main Injector) high intensity neutrino beam is generated at Fermilab by ex- 
tracting 120 GeV protons from the Main Injector and striking a graphite target. This produces pions and 
kaons which are focused into a 675 m long decay pipe by two magnetic focusing horns. Mesons travel down the 
decay pipe, filled with Helium at 0.9 atm. Undecayed hadrons and muons coming from the decay of mesons are 
absorbed and monitored. The resulting beam is made purely of neutrinos and antineutrinos. 

MINOS measures neutrino and antineutrino oscillations directly by utilizing NuMFs ability produce a v^- 
or F^-enhanced beam. In neutrino-mode the magnetic field in the horns is adjusted so that positively charged 
mesons are focused into the decay pipe which decay into neutrinos. In antineutrino-mode negative mesons are 
selected which enhances the antineutrino composition of the beam, as well as resulting in a pronounced peak 
at the energy range where oscillations are expected at around 3GeV. 



The analysis is performed in several steps: Selection of charged current neutrino and antineutrino events, 
measurement of the Near Detector and Far Detector energy spectra, calculation of the predicted Far Detector 
spectrum from the near, and a combined fit to neutrino and antineutrino datasets to measure common oscillation 
parameters and the NSI parameter, e. 



A preselection is applied to all data samples to ensure data quality. The same set of data quality criteria 
are required for both neutrino and antineutrino samples. We select on good beam type, good horn current and 
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good detector coil current as the running conditions require. We require that events have at least one track and 
have track vertices in the fiducial volume. 

Once data quality is satisfied, there are two major sources of background: neutral current events and events 
of the opposite charge sign. To separate charged current events from neutral currents, we employ a fc-Nearcst 
Neighbors (fcNN) algorithm. It is a multivariate method in which an N-dimensional space of track variables 
is constructed. A query event is identified as signal-like or background-like by its fc neighbors' characteristics. 
The distance in fcNN space is determined by 



D = 



\X, MC -Xf 



(9) 



where MC denotes the Monte Carlo events of known type (charged or neutral current) and Q is the query 
event. The particle identification parameter (fcNN PID) is then given by 



fcNN m = 



ks 



(10) 



where S and B denote signal and background respectively. In the cartoon in figure [3J red dots denote signal 
and blue dots background. In this case fc = 12, ks — 8, — 4: which gives the query event marked by the star 
a fcNNio of f. 
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FIG. 2: A cartoon demonstrating the use of the fcNN particle identification in two-dimensional fcNN space. 12 nearest 
neighbors lie in the dashed circle, red representing signal and blue representing background. This gives the event a fcNN 



PID value of § 



To select antineutrinos, we keep events with fcNN > 0.3. The fcNN PID distribution of antineutrinos (before 
selection) is shown in figure [3l We further require that events are reconstructed with positive charge. The 
efficiency and purity of this selector are shown in figure |4] 

The neutrino selection also uses the fcNN PID to separate charged current and neutral current events. An 
additional fcNN selector which improves low energy selection efficiency is also applied to select neutrino events. 
The overall neutrino CC/NC separator is a logical OR of these two discriminators. Events are required to have 
negative reconstructed charge. The efficiency and purity of the neutrino selection are shown in figure [5j 



Proceedings of the DPF-2011 Conference, Providence, RI, August 8-13, 2011 



5 




FIG. 3: The feNN PID distribution of antincutrino events in the Near Detector. Events above 0.3 are kept. 
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FIG. 4: Performance of the antineutrino selection (kNNio > 0.3, qp > 0). The dashed lines show the contamination 
before selection and the solid show efficiency and contamination after selection. 



B. Extrapolation 



The Near and Far Detector spectra are similar but not identical even in the absence of oscillations. Due to 
its proximity to the decay pipe, the Near Detector accepts boosted events that decay further downstream in 
the pipe whereas the Far Detector views it as a point source. To obtain an accurate prediction which takes into 
account these effects, a beam matrix is used. This matrix, shown in figure |6l is obtained by tracing a parent 
pion of known energy from Monte Carlo to its daughter neutrinos in the Near and Far Detectors. This matrix 
then muultiplies the measured Near Detector data spectrum to obtain a Far Detector prediction in a data- 
driven manner. Neutrinos and antineutrinos are extrapolated separately with their respective beam matrices, 
the figure shows the antineutrino matrix. 



C. Far Detector Data and Fitting 



We expect 2073 neutrino and 156 antineutrino events at the Far Detector in the absence of oscillations, and 
observe 1654 and 97 respectively. We fit for non-standard interactions in 100 bins per sample (for a total of 400 
bins in three neutrino and one antineutrino run periods). The effect of the four largest systematic uncertainties 
are included in the fit as penalty terms. Figure [7] shows the Far Detector energy spectra of neutrinos and 
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FIG. 5: Performance of the neutrino selection. The dashed lines show the contamination before selection and the solid 
show efficiency and contamination after selection. 
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FIG. 6: The beam matrix uses Monte Carlo truth information to trace a parent pion to its daughter neutrinos in the 
Near and Far Detectors. The matrix is then applied to the observed ND data spectrum to obtain the FD predicted 
spectrum. 



antincutrinos. 

The fit to non-standard interactions yields Am 2 = 2.56±0.15 x 1CT 3 eV 2 , sin 2 29> 0.90, and e=-0.187±0.16. 
The confidence intervals are shown in figure 
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FIG. 7: Far Detector reconstructed neutrino energy spectra of neutrino (left) and antineutrino samples. The red line 
shows the predicted spectrum without oscillations while the blue and green curves are fits to non-standard interactions 
and vacuum oscilations respectively. The bottom plots show the ratio of oscillated to unoscillated spectra. 
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FIG. 8: Allowed regions for Am 2 , sin 2 28, and e. Contours assume Gaussian statistics and include the effect of systematics. 
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